Frequency shift of cesium clock transition due to blackbody radiation 



E. J. Angstmann, V. A. Dzuba, and V. V. Flambaum 
School of Physics, University of New South Wales, Sydney 2052, Australia 
(Dated: February 2, 2008) 

We have performed ab initio calculations of the frequency shift induced by a static electric field 
on the cesium clock hyperfine transition. The calculations are used to find the frequency shifts due 



to blackbody radiation. Our result {5v/E 



-2.26(2) X 10"^"Hz/(V/m)^) is in good agreement 



with early measurements and ah initio calculations performed in other groups. We present argu- 
ments against recent claims that the actual value of the effect might be smaller. The difference (~ 
10%) between ab initio and semiempirical calculations is due to the contribution of the continuum 
spectrum to the sum over intermediate states. 

PACS numbers: 32.60.-(-i,31.30.Gs,31.25.Eb 



Atomic clocks are now important for both practical 
applications and fundamental physics. One of the domi- 
nant uncertainties in high-precision measurements of fre- 
quencies in atomic clocks is the ac stark shift induced 
by blackbody radiation (see e.g. There is some 

disagreement on the value of this shift. Early measure- 
ments 0, 0, Q and ah initio calculations 0, @] support 
a value which is close to —2.2 x 10~-^°Hz/(V/m)^ while 
more recent measurements [Tollllj | and semiempirical cal- 
culations 3, 9, 12] claim that actual number might be 
about 10% smaller. 

In the present work we have performed fully ah ini- 
tio calculations of the radiation frequency shift and have 
identified the source of the disagreement between differ- 
ent theoretical results as the contribution of the contin- 
uum spectrum states into summation over the complete 
set of intermediate states. The continuum spectrum was 
included in all the ah initio calculations and missed in 
the semiempirical considerations. We demonstrate that 
adding the contribution of the continuum spectrum to 
where it was missed brings all theoretical results to good 
agreement with each other and with early measurements. 

Blackbody radiation creates a temperature dependent 
electric field, described by the Planck radiation law 
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TT exp(ci;/fcT) — 1 



(1) 



This leads to the following expression for the average 
electric field radiated by a black body at temperature T: 



{E"^) = (831.9V/m)2[T(K)/300]^ 



(2) 



This electric field causes a temperature-dependent fre- 
quency shift of the atomic microwave clock transitions. 
It can be presented in the form (see, e.g. |]J) 



5u/uo = P{T/Tof [1 + e{T/Tof 



(3) 



Here Tq is usually assumed to be room temperature (Tq 
300iC). The frequency shift in a static electric field is 



5v ^ kE^ 



(4) 



Coefficients k and /3 are related by 



13 = — (831.9V/m)^ 

= fc X 7.529 X 10"^(V/m)2Hz"^ (for Cs), 



(5) 



while e is a small correction due to frequency distribution 
(^. In the present work we calculate the coefficient k. 

In the case when there is no other external electric field 
the radiation shift can be expressed in terms of the scalar 
hyperfine polarizability of the atom. This corresponds 
to averaging over all possible directions of the electric 
field. The hyperfine polarizability is the difference in the 
atomic polarizabilities between different hyperfine struc- 
ture states of the atom. The lowest-order effect is linear 
in the hyperfine interaction and quadratic in the electric 
field. The corresponding third-order perturbation theory 
expressions, after angular reduction have the form 
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m,j ^ ™Pj ' 
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Here 



= ^(2F' + 1)[F'(F' + 1) - /(/ + 1) - j{j + 1)] 

F' 

x|l/2 F'^\I-J\,I + J, 
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Ans is the hfs constant of the ns state, Am,n is the off- 
diagonal hfs matrix element, I is nuclear spin, F = I + J, 
J is total electron momentum of the atom in the ground 
state (J — 1/2), and j is total momentum of virtual p- 
states {j — 1/2,3/2). Summation goes over a complete 
set of ns, mpi/2 and states. 

In order to calculate frequency shift to the hfs transi- 
tions due to the electric field one needs to have a com- 
plete set of states and to have the energies, electric dipole 
transition amplitudes and hyperfine structure matrix el- 
ements corresponding to these states. It is possible to 
consider summation over the physical states and to then 
use experimental data to perform the calculations. The 
lowest valence states for which experimental data is usu- 
ally available dominate in the summation. Off-diagonal 
hfs matrix elements can be obtained to a high accuracy 
as the square root of the product of corresponding hfs 
constants: Am,n = VAZA^ (see, e.g. 0). However, 
the accuracy of this approach is limited by the need to 
include the tail contribution from highly excited states 
including states in the continuum. This contribution can 
be very significant and its calculation is not easier than 
the calculation of the whole sum. 

Therefore, in the present work we use an ab initio ap- 
proach in which high accuracy is achieved by the inclu- 
sion of all important many-body and relativistic effects. 
We make only one exception toward the semiempirical 
approach. The frequency shift is dominated by the term 
© which is proportional to the hfs in the ground state. 
It is natural to use experimental hfs constant in the domi- 
nating term to have more accurate results. Note however 
that the difference with complete ab initio calculations is 
small. 

Calculations start from the relativistic Hartree-Fock 
(RHF) method in the V^^^ approximation. This means 
that the initial RHF procedure is done for a closed-shell 
atomic core with the valence electron removed. After 
that, states of the external electron are calculated in the 
field of the frozen core. Correlations are included by 
means of the correlation potential method 16|. We use 
two different approximations for the correlation poten- 
tial, E. First, we calculate it in the lowest, second-order 
of the many-body perturbation theory (MBPT). We use 
notation E'^^^ for the corresponding correlation potential. 
Then we also include into E two classes of the higher- 
order terms: screening of the Coulomb interaction and 
hole-particle interaction (see, e.g. [I?] for details). These 
two effects are included in all orders of the MBPT and 
the corresponding correlation potential is named 

To calculate E^'^^ we need a complete set of single- 
electron orbitals. We use the B-spline technique 0,13 
to construct the basis. The orbitals are built as linear 
combinations of 50 B-splines in a cavity of radius AOob- 
The coefficients are chosen from the condition that the 
orbitals are eigenstates of the RHF Hamiltonian Hq of the 
closed-shell core. The E(°°) operator is calculated with 
the technique which combines solving equations for the 
Green functions (for the direct diagram) with the summa- 



tion over complete set of states (exchange diagram) [T^ . 

The correlation potential E is then used to build a new 
set of single-electron states, the so-called Brueckner or- 
bitals. This set is to be used in the summation in equa- 
tions ©, O and ©. Here again we use the B-spline 
technique to build the basis. The procedure is very sim- 
ilar to the construction of the RHF B-spline basis. The 
only difference is that new orbitals arc now the eigen- 
states of the Hq + T, Hamiltonian, where E is either E*^^^ 
orE(°°). 

Brueckner orbitals which correspond to the lowest va- 
lence states are good approximations to the real physical 
states. Their quality can be checked by comparing exper- 
imental and theoretical energies. Moreover, their quality 
can be further improved by rescaling the correlation po- 
tential E to fit experimental energies exactly. We do this 
by replacing the Hq + E with the Hq + AE Hamiltonian in 
which the rescaling parameter A is chosen for each partial 
wave to fit the energy of the first valence state. The val- 
ues of A are A^ = 0.8 and Ap = 0.85 for E^^) and A^ = 0.99 
and Ap = 0.95 for E(°°). Note that the values are very 
close to unity. This means that even without rescaling 
the accuracy is very good and only a small adjustment to 
the value of E is needed. Note also that since the rescal- 
ing procedure affects both energies and wave functions, it 
usually leads to improved values of the matrix elements 
of external fields. In fact, this is a semiempirical method 
to include omitted higher-order correlation corrections. 

Matrix elements of the hfs and electric dipole opera- 
tors are found by means of the time-dependent Hartree- 
Fock (TDHF) method [Ullll. This method is equivalent 
to the well-known random-phase approximation (RPA). 
In the TDHF method, single-electron wave functions are 
presented in the form ~ ipQ + dip, where ipo is unper- 
turbed wave function. It is an eigenstate of the RHF 
Hamiltonian Hq'- {Hq — eo)'0o = 0- '^V' is the correc- 
tion due to external field. It can be found be solving the 
TDHF equation 

iHo~eo)d4'^-Se4'o-Ho-SV''-'4'o, (9) 

where Se is the correction to the energy due to external 
field {6e = for the electric dipole operator) , F is the op- 
erator of the external field {Hhfs or eE-r), and 5V^^^ is 
the correction to the self-consistent potential of the core 
due to external field. The TDHF equations are solved 
self-consistently for all states in the core. Then matrix 
elements between any (core or valence) states n and m 
are given by 

i^nlF + SV^'-^lij^). (10) 

The best results are achieved when tpn and tpm are 
Brueckner orbitals calculated with rescaled correlation 
potential E. 

We use equation H10|l for all hfs and electric dipole 
matrix elements in © , {Tj) , and © except for the ground 
state hfs matrix element in ijSJ) where we use experimental 
data. 
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To check the accuracy of the calculations we perform 
calculations of the hfs in the ground state and of the 
static scalar polarizability. Polarizability is given by the 
expression 

2^ |(a|H|m)p 
ao(a) = -2^ Z (11) 

m ^™ 

which is very similar to the term ||HJ| for the frequency 
shift. The most important difference is that the energy 
denominator is squared in term (|SJ but not in . This 
means better convergence with respect to the summa- 
tion over complete set of states for term (jHl than for 
Therefore, if good accuracy is achieved for polariz- 
abilities, even better accuracy should be expected for the 
term © (see also Ref. ^). 

However, the behavior of the other two terms, © and 
O, is very different and calculation of polarizabilities 
tells us little about accuracy for these terms. Therefore, 
we also perform detailed calculations of the hfs constants 
of the ground state. Inclusion of core polarization (second 
term in (|10|l ') involves summation over the complete set of 
states similar to what is needed for term ^ . Comparing 
experimental and theoretical hfs is a good test for the 
accuracy for this term. 

In addition to term l|10(l . we also include two smaller 
contributions to the hfs: structure radiation and the cor- 
rection due to the change of the normalization of the 
wave function. Our final result for the hfs constant is 
2278 MHz which is in excellent agreement with the ex- 
perimental value 2298 MHz |2^. The result for static po- 
larizability is q;o = 399.0 Oq which is also in a v ery good 
agreement with experimental value 401.0(6) [lj|- 

Table m presents contributions of terms ©, Q and ^ 
into the total frequency shift of the hfs transition for the 
ground state of ^■^^Cs calculated in different approxima- 
tions. Term ^ dominates while term (Q is small but 
still important. Results obtained with E^^' and 
differ significantly (14%). However, after rescaling the 
results for both E^^^ and E(°°) come within a fraction of 
a per cent of each other. Naturally, rescaling has a larger 
effect on results obtained with fP^y This means that the 
rescaling really imitates the effect of higher-order corre- 
lations and should lead to more accurate results. 

In summary, we have three ways of estimation of the 
accuracy of calculations: (a) calculation of static polariz- 
ability (0.5% accuracy); (b) calculation of the hfs (0.9% 
accuracy); and (c) comparision of the results obtained in 
different most accurate approximations (three last lines 
of Table nj, which differ by about 0.3%. Therefore, we 
can say that the accuracy of the calculations is not worse 
that 1%. Our final result is 

k = -2.26(2) X 10"^°Hz/(V/m)2. (12) 

This corresponds to /3 = -1.70(2) x lO"^"'. To obtain 
frequency shift at finite temperature one needs to sub- 
stitute this value into equation ||3J|. For accurate results 
one also needs to know the value of e. It was estimated in 



TABLE I: Contribution of terms ®, 0, and © to the fre- 
quencies of the hyperfine transition in the ground state of 
^^'^Cs {Svo/E'^ X 10"^° Hz/(V/m)2) in different approxima- 
tions. 



E 




Q 


© 


Total 


f;(2)a 


-0.9419 


0.0210 


-1.0722 


-1.9931 




-1.0239 


0.0229 


-1.2688 


-2.2697 


^ (oc) c 


-1.0148 


0.0232 


-1.2706 


-2.2622 




-1.0167 


0.0230 


-1.2695 


-2.2632 



"S^^' is the second-order correlation potential. 
''Rescaled t. 

<^Yi(°°) is the all-order correlation potential. 



Ref. m in single-resonance approximation and found to 
be 0.014. In many-resonance calculation e will be 10-20% 
smaller. 

We present our final result for the frequency shift to- 
gether with other theoretical and experimental results in 
Table IhI Our value is in good agreement with early mea- 
surements S 3.nd ah initio calculations 0, Q while 
recent measurements [Tol lllj and semiempirical calcula- 
tions 007 give the value which is about 10% smaller. 
Less accurate measurements of Bauch and Schroder 
cover both cases. We cannot comment on disagreement 
between experimental results. However, the source of dis- 
agreement between theoretical results seems to be clear. 
It comes from the contribution of the continuum spec- 
trum to the summation over the complete set of states 
in term ©. This term has off-diagonal hfs matrix el- 
ements between the ground state and excited states. 
Since the hfs interaction is localized over short distances 
(~ ao/Z) it emphasizes the contribution of states with 
high energies including states in the continuum (since 
ApAx > ?i, a small area of localization (Aa:) allows high 
momentum (p) and thus high energy). In our calcula- 
tions the contribution of states above 7p into term (jSJ is 
-0.35x 10-iHz/(V/m)2 which is 15% of the total answer. 

In contrast, states above 7p contribute only about 
0.05% of the total value of term JHJ. This is because 
the summation goes over the matrix elements of the elec- 
tric dipole operator which is large on large distances and 
thus suppresses the contribution of high-energy states. It 
is not surprising therefore that a semiempirical consider- 
ation, involving only discrete spectrum states, gives very 
good results for the atomic polarizabilities (see, e.g. 0). 
However, let us stress once more that the calculation of 
polarizabilities checks only term and tells us very lit- 
tle about the accuracy of other two terms, ® and Q). 

The contribution of the states above 7p is even more 
important for term ((7|). Their contribution is about 30% 
of the total value of this term. However, the term itself 
is small and its accurate treatment is less important. 

In ah initio calculations by Lee et al 01 summation over 
complete set of states is reduced to solving of a radial 
equation (equations of this type are often called Stern- 
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TABLE II: Electrostatic frequency shifts for the hyperfine 
transition of Cs {Svo/E'^ x 10~^° Hz/(V/m)^) ; comparison 
with other calculations and measurements. 



This 


Other 


Ref 


Measurements 


Ref 


work 


calculations 








2.26(2) 


-1.9(2) 




-2.29(7) 


[2] 




-2.2302 




-2.25(5) 






-2.28 




-2.17(26) 






-1.97(9) 


m 


-2.271(4) 


m 




-2.06(1) 


[12j 


-1.89(12) 


m 




-2.268(8) 


113J 


-2.03(4) 





heimer equation after one of the authors of this work). 
This approach does include the contribution of the con- 
tinuum spectrum and the result is in very good agreement 
with ours (see Table HDl. 

In other ah initio calculations by Pal'chikov et al Q 
summation is done via Green functions. This corresponds 



to summation over the complete set of states and does 
include the continuum spectrum. Again, the result is in 
very good agreement with other ah initio calculations (0 
and the present work). 

Recent calculations by Beloy et al applied a mixed 
approach, with extensive use of experimental data for 
lower cesium states and ah initio summation over higher 
states including continuum. The result is in good agree- 
ment with fully ah initio calculations. 

In contrast, analysis performed in 0,0, ^| is limited 
to discrete spectrum. Adding —0.34 x 10~^Hz/(V/m)^ 
(which is total tail contribution from all three terms 
(0 and ijHI) found in our calculation) to the results of 
Feitchner et al Q and Micalizio et al Q brings them to 
excellent agreement with ah initio calculations. The same 
modification of the result by Ulzega et al makes it 
a little bit too large but still closer to other results than 
without the tail contribution. 

We are grateful to S. Ulzega, W. Itano, and A. Dere- 
vianko for useful comments and references. 
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